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Abstract 

This paper presents measurements of K° and A production in neutral current, deep inelastic 
scattering of 26.7 GeV electrons and 820 GeV protons in the kinematic range 10 < Q 2 < 640 
GeV 2 , 0.0003 < x < 0.01, and y > 0.04. Average multiplicities for K° and A production are 
determined for transverse momenta pt > 0.5 GeV and pseudorapidities |ry| < 1.3. The multiplicities 
favour a stronger strange to light quark suppression in the fragmentation chain than found in e + e~ 
experiments. The production properties of K°'s in events with and without a large rapidity gap 
with respect to the proton direction are compared. The ratio of neutral K 0, s to charged particles per 
event in the measured kinematic range is, within the present statistics, the same in both samples. 
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1 Introduction 



The investigation of strange particle production in neutral current, deep inelastic scattering (DIS) 
interactions could provide information about the s-quarks in the nucleon, about the boson-gluon 
fusion process and, above all, the parton fragmentation process. Strange particle production has been 
measured previously by experiments where the 7*p centre-of-mass energy, W, is at least one order of 
magnitude lower than at HERA jl], ^, [|, [|] . The ratio of strange particle to light non-strange particle 
production of approximately 1:5 is ascribed to a reduced probability of strange quark creation in the 
parton fragmentation chain. In simulation programs based on the Lund scheme it is parametrised by 
the strange quark suppression factor P s /P u . Here P s and P u are the probabilities for creating s— or 
u, d— quarks from the vacuum during the fragmentation process. A detailed review of our knowledge 
on heavy quark suppression is given in Q. In hadron-hadron collisions an increasing P s /P u is found 
with increasing centre-of-mass energy. Also indications of a dependence of the strangeness suppression 
factor on the region of phase space under investigation are reported. The values found vary between 
about 0.15 and 0.55 with a mean value close to 0.3 (see for example [^, 0). The parameters for the 
hadronisation process in the present day electron-proton Monte Carlo event generators are obtained 
from fits to e + e~ data and are assumed to be the same in DIS experiments due to jet universality. 
The longitudinal phase space of the 7*p interactions at HERA can be divided into three regions where 
different processes are expected to dominate. These processes also appear in e + e~ or hadron-hadron 
scattering: a) the fragmentation region of the struck quark, which resembles that of one of the pair- 
produced quarks in e + e~ annihilation experiments; b) the fragmentation of the proton remnant, which 
resembles the fragmentation in hadron colliders; and c) the hadronic centre-of-mass central rapidity 
region, where the colour flow between the struck quark and the proton remnant evolves. The latter 
region exists in both e + e~ and hadron collider experiments. The acceptance of our central tracking 
detector allows us to study K° production in the fragmentation region of the struck quark and the 
central rapidity region. The part of the event which is well inside our detector acceptance is dominated 
by particles originating from the central rapidity region. 

In about 10% of the DIS events no proton remnant is detected in the ZEUS detector, resulting in a 
large rapidity gap (LRG) between the acceptance limit in the proton direction and the first visible 
particle in the detector [||, |9|. The properties of these events are consistent with the assumption 
that the exchanged photon is scattered off a colourless object emitted by the proton. This object is 
generically called a pomeron. There exist indications that the pomeron has a partonic substructure 
|To| , |8| but the nature of its constituents is still under investigation. A natural assumption is that they 
are quarks and gluons or a combination of both. It is expected that the strange quark content of the 
pomeron could affect the strange particle multiplicity in the final state of these events. 
The investigation of strange particle production allows us to connect results from e + e~ experiments 
and from hadron collider experiments. This paper is a first step of such a program. We compare 
the K° and A multiplicities^] and their momentum and angular distributions in the new kinematic 
region of HERA with extrapolations from Monte Carlo models based on the results of lower energy 
experiments. The Q 2 evolution of the K° multiplicity is studied. The production of K 0, s in events 
with a large rapidity gap is compared to that of events without a large rapidity gap. 
All studies are performed in the HERA laboratory frame and are restricted to a kinematic range where 
the tracking acceptance is high and well understood. 



1 Throughout this paper, a reference to a particle includes a reference to its antiparticle. 
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2 Experimental setup 



HERA machine conditions 

The data were collected at the electron-proton collider HERA using the ZEUS detector during the 
1993 running period. HERA collided 26.7 GeV electrons with 820 GeV protons. 84 bunches were filled 
for each beam and in addition 10 electron and 6 proton bunches were left unpaired for background 
studies. The typical electron and proton currents were 10 mA leading to a typical instantaneous 
luminosity of 6 • 10 29 cm~ 2 s . An integrated luminosity of 0.55 pb _1 of data was collected in 1993. 

The ZEUS detector 



ZEUS is a multipurpose, magnetic detector which has been described elsewhere [11]. Here we give a 
brief description concentrating on those parts of the detector relevant for the present analysis. 
Charged particles are tracked by the inner tracking detectors which operate in a magnetic field of 1.43 
T provided by a thin superconducting solenoid surrounding the tracking detectors. Immediately after 
the beampipe there is a cylindrical drift chamber, the vertex detector (VXD), which consists of 120 



radial cells, each with 12 sense wires |12]. The achieved resolution is 50 /im in the central region of a 
cell and 150 fim near the edges. Surrounding the VXD is the central tracking detector (CTD) which 
consists of 72 cylindrical drift chamber layers, organised into 9 "superlayers" |13|]. Each superlayer 
consists either of wires parallel (axial) to the beam axis or of wires inclined at a small angle to give a 
stereo view. With the present understanding of the chamber, a spatial resolution of 260 /jm has been 
achieved. The hit efficiency of the chamber is greater than 95%. 

In events with charged particle tracks, using the combined data from both chambers, reconstructed pri- 
mary vertex position resolutions of 0.6 cm in the Z direction and 0.1 cm in the XY plane are measuredpl. 



The resolution in transverse momentum for full length tracks is ct(pt)/pt = V (0.005 pt) 2 + (0.016) 2 
(pr in GeV). 

The solenoid is surrounded by a high resolution uranium-scintillator calorimeter divided into three 
parts, forward (FCAL), barrel (BCAL) and rear (RCAL). Holes of 20 x 20 cm 2 in the centre of FCAL 
and RCAL are required to accommodate the HERA beam pipe. Each of the calorimeter parts is 
subdivided into towers which in turn are segmented longitudinally into electromagnetic (EMC) and 
hadronic (HAC) sections. A section of a tower is called a cell and is read out by two photomultiplier 



tubes. A detailed description of the calorimeter is given in [14|. 

For measuring the luminosity via the Bethe-Heitler process ep — > e'p'^, as well as for tagging very 
small Q 2 processes, two lead-scintillator calorimeters are used [15|. Bremsstrahlung photons emerging 



from the electron-proton interaction region at angles 6'^ < 0.5 mrad with respect to the electron beam 
axis hit the photon calorimeter at 107 m from the interaction point (IP). Electrons emitted from the 
IP at scattering angles less than or equal to 6 mrad and with energies between 20% and 90% of the 
incident electron energy are deflected by beam magnets and hit the electron calorimeter placed 35 m 
from the IP. 



3 HERA kinematics 

The kinematics of deep inelastic scattering processes at HERA, e~p — > e~h, where h is the hadronic 
final state, can be described by the Lorentz invariant variables Q 2 , x and y. Here — Q 2 is the square 
of the four-momentum transfer between the incoming electron and the scattered electron; x, in the 
naive quark-parton model, is the fractional momentum of the struck quark in the proton, and y is 
the relative energy transfer of the electron to the hadronic system. The variables are related by 
Q 2 = sxy, where s is the squared invariant mass of the ep system. Q 2 , x and y can be calculated 

2 The ZEUS coordinate system is denned as right-handed with the Z axis pointing in the proton beam direction and 
the X axis horizontal pointing towards the centre of HERA. The polar angle 6 is defined with respect to the Z-direction. 
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from the kinematic variables of the scattered electron, from the hadronic final state variables, or from 
a combination of both. The optimal reconstruction method depends on the event kinematics and the 
detector resolution. 

In this paper we use the double angle method [|^] to calculate the Q 2 and x variables: 

n 2 Aw i sin 7h (l + cos(9 e ) 



XDA 



sin j h + sin 9 e - sm.{~f h + 6 e ) ' 
E e sin7 fe + sin6* e + sin (7^ + 9 e ) 
E p sin 7/j + sin 9 e - sin(7/> + 9 e )' 



Here E e and £" p are the initial electron and proton energies; 9 e is the electron scattering angle with 
respect to the incident proton direction and is the polar angle of a massless object balancing the 
momentum vector of the scattered electron to satisfy four-momentum conservation. In the naive 
quark-parton model 7/, is the scattering angle of the struck quark. It is determined from the hadronic 
energy flow in the calorimeter: 

rn _ (Epx)1 + (Epy)1-(EE-pz)1 
7/1 (EPx)l + (Epy)l + (EE- P z) 2 h 

Here the sums run over all calorimeter cells which are not assigned to the scattered electron and 
(px,Py,Pz) is the momentum vector assigned to each cell of energy E. The cell angles are calculated 
from the geometric centres of the cells and the vertex position of the event. 

Using the hadronic energy flow of the final state, y can be calculated according to the Jacquet-Blondel 
method [17|: 

1 

VJB = Yl (E-Pz)h- 

2Ee h 

For background rejection we also calculate y using the electron information: 

E' 

Ve = l ~2E~ <yl ~ cos9e ^ 

where E' e is the energy of the scattered electron. The square of the centre-of-mass energy of the virtual 
photon-proton system, 7*p, is calculated using: 

Wh A = m 2 p + Q 2 DA (—-l), 

where m p is the proton mass. We use the described methods for calculating the kinematic variables 
and do not mention them explicitly below except when necessary. 



4 Event selection 
4.1 Trigger conditions 

The trigger is organised in three levels |fLl] ]. For DIS events, the first level trigger (FLT) requires at 
least one of three conditions of energy sums in the EMC calorimeter: the BCAL EMC energy exceeds 
3.4 GeV; the RCAL EMC energy (excluding the innermost towers surrounding the beam pipe) exceeds 
2.0 GeV; or the RCAL EMC energy (including those towers) exceeds 3.75 GeV. 

The second level trigger (SLT) rejects proton beam-gas events by using the times measured in the 
calorimeter cells. The DIS trigger rate of the SLT is about one-tenth of the FLT DIS trigger rate. 
The loss of DIS events at the SLT is negligible. 
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The third level trigger (TLT) has the full event information available and uses physics-based filters. 
It applies tighter timing cuts to suppress beam-gas background further and also rejects beam halo 
muons and cosmic muons. The TLT selects DIS event candidates by calculating: 

5= J2 Ei ' > 20 GeV - 2 E 7 , 

i 

where E^ and 6i are the energy and the polar angle of the energy deposits in the central calorimeter. 

is the energy measured in the photon calorimeter of the luminosity monitor. The summation runs 
over all energy deposits in the calorimeter cells. For fully contained DIS events 5 ~ 2E e = 53.4 GeV. 
Photoproduction events have low values of 5 compared to DIS events because the scattered electron 
escapes in the beam pipe hole of the rear calorimeter. For events with Q 2 less than ~ 4 GeV 2 the 
calorimeter cannot detect the scattered electron. 

For events with the scattered electron detected in the calorimeter, the trigger acceptance was essen- 
tially independent of the DIS hadronic final state. It was greater than 97% for Q 2 > 10 GeV 2 and 
independent of Q 2 . A total of 7 • 10 6 events passed the TLT and was written to tape during the 1993 
running period. 



4.2 Offline event selection 



The offline selection of DIS events is similar to that described in our earlier publication [18]. The 
characteristic signature of a DIS event is the scattered electron detected in the uranium scintillator 
calorimeter. The pattern of energy deposition in the calorimeter cells is used to identify an electron 
candidate. We use the following criteria to select a sample of DIS events: 

• a scattered electron candidate has to be found with E' e > 5 GeV and an impact point at the 
RCAL surface outside a square of 32 x 32 cm 2 centred on the beam line. This requirement 
ensures that the electromagnetic shower is fully contained within the calorimeter and its impact 
point can be reconstructed with sufficient accuracy; 

• y e < 0.95 to reduce photoproduction background; 

• 35 GeV < 5 < 60 GeV to remove photoproduction events and to suppress events with hard 
initial state radiation; 

• -50 cm < Z < 40 cm, where Z is the position of the event vertex reconstructed from the CTD. 
This requirement rejects beam-gas and cosmic ray events. 

From Monte Carlo studies we find an average electron finding efficiency of 95% in the kinematic range 
considered, being above 98% in most of the kinematic range and dropping below 70% for high y events. 
The purity is better than 96% for electron energies above 10 GeV and drops to about 60% at high y. 
A total of 91000 events survive these criteria. 

The particle multiplicity and the kinematics of particle production depend on Q 2 and x. We have 
restricted our analysis to a kinematic range in Q 2 , x and y in which migration effects are small (|18|| 
and have little influence on the momentum and angular distributions of the K 0, s and A's. We chose 
the following range: 

• 10 GeV 2 < Q 2 < 640 GeV 2 ; 

• 0.0003 < x < 0.01; 

• y > 0.04. 

The Q 2 and x variables are calculated according to the double angle method and y with the Jacquet- 
Blondel method. After applying these criteria to the previously selected sample we are left with 27500 
events. 



4 



5 Monte Carlo simulation 



Monte Carlo event simulation is used to determine the acceptance and resolution of the ZEUS detector. 
The simulation is based on the GEANT 3.13 |l9| program and incorporates the knowledge of the 
detector and the trigger. 

Simulation of normal DIS events 

Neutral current DIS events with Q 2 > 4 GeV 2 were generated using the HERACLES 4.4 program 
which incorporates first order electroweak corrections. The Monte Carlo program LEPTO 6.1 



interfaced to HERACLES via the program DJANGO 6.0 [22], was used to simulate QCD cascades 



and fragmentation. The parton cascade was modelled in different ways: 

• the colour-dipole model including the boson-gluon fusion process (CDM) as implemented in the 



ARIADNE 4.03 [23| program was used. In this model coherence effects are implicitly included 



in the formalism of the parton cascade; 

• matrix element calculations plus the parton shower option (MEPS) as implemented in LEPTO 
were used, where coherence effects in the final state cascade are included by angular ordering of 
successive parton emissions. 

These models use the Lund string fragmentation [24| for the hadronisation phase as implemented in 
JETSET 7.3 || . 



For the CDM event sample the MRSD'_ parton density parametrisation for the proton was used [26|. 



The GRV [27] parametrisation was used for the MEPS data set. These parametrisations describe 



reasonably the HERA measurements of the proton structure function F2 [28, 29] 



The simulations predict that about 10% of the Kg's are produced in charm events and about 5% 
originate from sea quarks in the proton. The remaining ~ 85% of the Kg's are created in the 
fragmentation chain depending on the actual value of the strange-quark suppression factor P s /P u . 
The parameters of the Monte Carlo models are set to their default values (P s /P u = 0.3). We have also 
generated events with P s /P u = 0.2 as suggested in || and we have compared the predictions of the 
simulations with the measured rates. Since the MEPS model and the CDM model behave similarly 
when reducing the P s /P u parameter, we only show the predictions with P s /P u = 0.2 for the CDM 
model. 

Simulation of large rapidity gap DIS events 



Our previous study M shows that diffractive models, specifically POMPYT [30] and a model by 



Nikolaev and Zakharov [31] as implemented in our Monte Carlo program NZ p^[ , give adequate 
descriptions of the properties of the LRG events. We have used POMPYT and NZ event samples 
for our study of K° multiplicities in events with a large rapidity gap. The POMPYT Monte Carlo 



program uses an implementation of the Ingelman and Schlein model [33], describing high energy 
diffractive processes. In this model the virtual photon interacts with the constituents of the pomeron 
emitted by the proton. Factorisation is assumed in the sense that the pomeron emission and the 
pomeron structure are independent. The current version of POMPYT contains no strange quark 
constituents for the pomeron. The NZ Monte Carlo model on the other hand is non-factorisable. Here 
the virtual photon fluctuates into a qq or a qqg state and interacts with a colourless two-gluon system 
emitted by the proton. The qqg states were fragmented as if they were qq states and the flavours are 
generated in 90% of the cases as (u, d) and in 10% as s. 
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6 Selection of Kg and A candidates 



K s particles are identified in the decay channel Kg — > tt + it~ and A particles are detected in the 
channel A — > p-7r~. Due to their lifetime of O(10 _10 s) and their typical momenta of about 1 GeV they 
have an average decay length of a few centimetres, which results in secondary vertices well separated 
from the primary event vertex. 

Tracks are reconstructed using the CTD and the VXD. The track finding algorithm starts with hits in 
the outermost axial superlayers of the CTD. As the trajectory is followed inwards to the beam axis, 
more hits on the axial wires and from the VXD are assigned to the track. The resulting circle in 
the transverse plane is used for the pattern recognition in the stereo superlayers. The momentum is 
determined in a 5-parameter helix fit. Multiple Coulomb scattering in the beam pipe and in the outer 
walls of the VXD is taken into account in the evaluation of the covariance matrix. 
The primary event vertex is determined from a \ 2 fit performed with the tracks using the perigee 
parametrisation [Q| and assuming that the tracks come from a common point in space. A track is 
considered not to be associated with the primary vertex if the \ 2 f° r the primary vertex fit increases 
significantly when the track is included in the fit. 

The systematic effects in the CTD are most serious for low tracks and for tracks which traverse the 
inhomogeneous part of the magnetic field at the ends of the CTD. The reconstructed tracks used in 
this analysis were required to have a transverse momentum pi > 0.2 GeV and a polar angle between 
25° < 6 < 155°. In terms of pseudorapidity, rj = — log (tan (0/2)), this corresponds to \rj\ < 1.5. This 
is the region where the CTD response and systematics are well understood. 

6.1 K s Identification 

To search for Kg, we examine pairs of oppositely charged tracks to find a secondary vertex. We refer 
to these tracks as daughter tracks. At least one of the daughter tracks is not allowed to be associated 
with the primary vertex and track pairs which do not intersect when projected into the transverse 
plane are rejected. 

For each remaining track pair, we obtain the momentum of the Kg candidate by calculating the 
momenta of the individual tracks at their intersection point and adding them. K s candidates with 
transverse momenta below 0.5 GeV or above 4 GeV or with directions of flight too near to the beam 
pipe, \q\ > 1.3, are removed. 

The background in the mass region of the Kg is reduced by applying the following criteria: 

• cos(axy) > 0.99, where axY is the angle in the transverse plane between the direction of flight 
of the Kg candidate and its reconstructed momentum direction; 

• the separation in Z between the two tracks at their XY intersection point has to be \AZ\ < 
2.5 cm. The coordinates of the Kg decay vertex are set to the XY coordinates of the intersection 
point of the track circles and the Z coordinate is chosen to be in the centre between the closest 
approaches in Z of the two track circles; 

• the proper lifetime of the candidates, cr = {LMc) /p, has to be less than 10 cm. Here L is the 
decay length, p is the momentum and M is the invariant mass of the candidate; 

• to reduce background arising from photon conversions into e + e~ pairs, pairs of tracks considered 
as electrons must have an effective mass M ee > 50 MeV (see Fig. [j]); 

• to eliminate A contamination of the Kg signal, candidates with a mass hypothesis M p7T < 
1.12 GeV are rejected (see Fig. |l|). 

Using these criteria (summarised in Tab. ||) we obtain the Kg signal shown in Fig. [2|a. We fit the 
7T + 7r - mass spectrum with a Gaussian and a linear background in the region 0.4 to 0.6 GeV. The 
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fitted mass is 497.4 ± 0.3 MeV and the standard deviation is 7.8 ± 0.3 MeV. The mass value and 
width of the signal are well reproduced by the Monte Carlo simulations. In the signal region we find a 
total of 971 Kg mesons on top of a background of about 150 7T7r-combinations. The K% signal region 
extends from 474 to 521 MeV. The average lifetime of the K s mesons was determined by fitting the 
exponential form exp(— ct/ct^q ) to the acceptance corrected cr lifetime distribution. Here the cr 

upper limit was relaxed to 20 cm and all other selection criteria were set to their default value. The 
result CTj^o = 2.66 ± 0.11 ± 0.06 cm is consistent with the world average of 2.676 ± 0.006 cm given 

5 

in |35|]. The systematic uncertainty includes the variation of the number of bins used in the fit and 
tightening or loosening the selection criteria. 



6.2 A Identification 

The A identification closely resembles the K$ identification. The daughter track with the higher 
momentum is considered to be the proton. No daughter track is allowed to be associated with the 
primary vertex. The cr upper limit is increased to 40 cm in order to account for the longer lifetime 
of the A. Requiring M nn < 0.481 GeV removes the background from Kg mesons (see Fig. [l]). Since 
there is no clear A signal seen for candidates with above 3.5 GeV, this value is chosen as the upper 
limit of the investigated momentum range. 

Figure §b shows the A signal obtained. We fit the pir mass spectrum from 1085 to 1185 MeV. The 
fit yields a mass of 1116.2 ± 0.4 MeV with a standard deviation of 3.0 ± 0.5 MeV. The Monte Carlo 
simulation reproduces well the A mass position and width. Within the signal region we find 80 A 
baryons and 18 background combinations. The signal region runs from 1107 to 1125 MeV. Of the 80 
A baryons, (60 ± 5)% are A and the remaining (40 ± 5)% are A. The determination of the average 
lifetime of the A from the lifetime distribution gives cr^ = 7.3 ± 2.2 ± 0.5 cm, consistent with the 



value of 7.89 cm given in [35] 



Selection parameters for candidates 


n s 


A 


cos(axy) 


> 0.99 


> 0.99 


\AZ\ [cm] 


< 2.5 


< 2.5 


cr [cm] 


< 10 


< 40 


M p7V [GeV] 


> 1.12 




M nw [GeV] 




< 0.481 


M ee [GeV] 


> 0.05 


> 0.05 


PT daughter— tracks [GeV] 


> 0.2 


> 0.2 


$ daughter— tracks [ ] 


[25, 155] 


[25, 155] 


No. of tracks from primary vertex 


< 1 





r\ range 


[-1.3, 1.3] 


[-1.3, 1.3] 


PT [GeV] range 


[0.5, 4.0] 


[0.5, 3.5] 



Table 1: Selection criteria for Ko and A identification. 



7 Data correction 

This analysis uses two types of selection criteria. The first kind is event based and selects a reasonably 
pure sample of DIS events with minimal contamination from background (photoproduction, beam-gas, 
cosmic-ray events). The second kind of selection criteria is particle based and selects a sample of K s 
and A particles from the event sample defined above. 
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We find a 90% event selection efficiency, where we define the efficiency as the ratio of the number of 
Monte Carlo events passing all the event selection criteria (including those that restrict the kinematic 
range in Q 2 ,x and y) to the total number of generated events in the restricted kinematic region. 
We have restricted the K° and A kinematic ranges to regions where our systematic uncertainties are 
small: their pseudorapidity is limited to — 1.3 < 77 < 1.3 and their transverse momentum is restricted 
to apr between 0.5 GeV and 4.0 GeV (3.5 GeV) for K 0, s (A's). We do not extrapolate our results to 
the full pt and rj range in order not to be dominated by model predictions. The models are known to 
have uncertainties especially in the low pi region and are not yet compared to particle properties in 
the proton fragmentation region of HERA events. 

The Kg and A reconstruction efficiencies were determined as a function of px and 77. For each particle 
type, the efficiency in a given (pr> v) Dm was defined as the ratio of the number of reconstructed 
particles in the bin to the number of generated particles in the bin. The rj and px resolutions are less 
than 5% of the bin width chosen for the plots and show no systematic shifts. The DIS Monte Carlo 
events that passed all the selection criteria were used for these calculations. The Kg reconstruction 
efficiency in the kinematic region considered varies between 20% for low px and 55% for px above 
1.5 GeV. The efficiency varies in 77 from 30% around 77 = ±1.3 to 40% for K s 's moving transversely 
to the beam direction (77 = 0). The A reconstruction efficiency varies between 5% for low transverse 
momentum and approaches 20% for high px- The efficiency varies in 77 between 10% and 15%. The 
largest loss of true KgS and A's results from the collinearity requirement (axy) and the requirement 
that daughter tracks are unassociated with the event vertex. Each requirement rejects about 25% of 
the candidates if no other selection criterion is applied. 

The if" (A) measurements are corrected for the above efficiencies as well as for the branching ratios 
K° to Kg and K s — > tt + tt~ (A — > pit) |^]. No corrections were made to the measurements for 
migrations and initial state radiation effects since the changes predicted from Monte Carlo studies are 
small. Instead we include these effects in our systematic error analysis (see section 9). 
The analysis procedure was checked using the reconstructed CDM (MEPS) Monte Carlo events as if 
they were data events and correcting them with the efficiencies obtained with the MEPS (CDM) sam- 
ples. The corrected Monte Carlo distributions agreed at the 5% level with the generated distributions. 
For the comparison of K° production in events with and without a large rapidity gap, the two- 
dimensional (pr, v) efficiencies were determined from the standard DIS Monte Carlo events satisfying 
the additional requirement W > 140 GeV (see section 8.2 for details). This corresponds to a restriction 
to y > 0.22. Both non-LRG (NRG) and LRG data samples were corrected with the same efficiencies. 
It has been checked that the corrected and generated K s distributions of the LRG Monte Carlo events 
agree well when using the efficiencies of those Monte Carlo sets. 

The ratio of if to charged particle multiplicity, N(if )/N (tracks), is investigated below. The charged 
particle multiplicity, N(tracks), is determined for charged particles originating at the primary vertex 
and produced in the restricted kinematic range [77! < 1.3 and px > 0.2 GeV. The number of recon- 
structed tracks is corrected for tracking inefficiencies, wrong assignments by the vertex finding routine 
to the decay products of long lived particles and pair conversions by using standard Monte Carlo 
techniques. The Monte Carlo corrections for particle based selection criteria were below 10%. 

8 Results 

8.1 K° and A multiplicity distributions 

Figure || shows the differential K° multiplicity as a function of px an d 77. The inner error bars are 
statistical errors and the outer ones statistical and systematic errors added in quadrature. The dis- 
tributions are normalised by the number of events N ev . The predictions of the CDM and the MEPS 
models are overlaid. The two curves for the CDM sample are generated with different strange-quark 
suppression factors P s /P u . The predicted multiplicity for the default strange-quark suppression factor 
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of 0.3 is higher than measured. Using the smaller suppression factor of 0.2 reduces the predicted mul- 
tiplicity to a value closer to that observed in the data. Both parameters give a reasonable description 
of the measured shapes. 

For events with 10 < Q 2 < 640 GeV 2 , 0.0003 < x < 0.01 and y > 0.04, the number of neutral 
kaons per event with 0.5 < pp < 4.0 GeV and \rj\ < 1.3 is 0.289 ± 0.015 ± 0.014. The first error 
is statistical, the second error is systematic. A function of the form Ci/pr • exp (C2pr) fits well the 
measured l/N ev ■ dN(K°)/dpj^ distribution as a function of pr over the pr range shown in Fig. [3|. 
C\ and C2 are constants. The slope, C2, of the pr distribution for the K°'s is —1.31 ± 0.09 ± 0.06 
GeV" 1 . These values, together with the predictions from Monte Carlo models, are listed in Tab. |2|. 
According to Monte Carlo studies, the fraction of K s produced in the restricted pr and r\ range is 
23% of the total number of K 0, s produced in the final state. 





N(K U ) / event 


pr slope [GeV : ] 


Data 

CDM 

with P s /P u = 0.3 
with P s /P u = 0.2 

MEPS 

with P s /P u = 0.3 


0.289 ± 0.015 ± 0.014 

0.342 ± 0.005 
0.264 ± 0.003 

0.348 ± 0.006 


-1.31 ± 0.09 ± 0.06 

-1.40 ± 0.05 
-1.37 ± 0.04 

-1.36 ± 0.05 



Table 2: Results of the K° measurement for events with 10 < Q 2 < 640 GeV 2 , 0.0003 < x < 0.01, 
y > 0.04 and for a K° with 0.5 < pr < 4.0 GeV and \r]\ < 1.3. The two CDM samples have been 
generated with a different strange-quark suppression factor P s /P u . 

Figures ^|a, b show the differential A multiplicity as a function of the transverse momentum and 
the pseudorapidity. The predictions of the CDM and the MEPS Monte Carlo are also displayed in 
Fig. |||. The two CDM curves correspond to samples generated with different strange-quark suppression 
factors P s /P u . The number of A's with 0.5 < pr < 3.5 GeV and |jy| < 1.3 produced per event is 
0.038 ± 0.006 ± 0.002 for events with 10 < Q 2 < 640 GeV 2 , 0.0003 < x < 0.01, y > 0.04. The 
measured slope of the pp distribution of the A is —1.4 ± 0.3 ± 0.1 GeV -1 , which, due to the large 
statistical uncertainty, is still in agreement with the model predictions. These values, together with 
the predictions of the models are listed in Tab. ||. Monte Carlo studies predict that 16% (25%) of the 
total number of A's (A's) will be inside this restricted pr and 77 region. 

The measured K° and A multiplicities seem to be better described by a model with a strangeness 
suppression factor of 0.2. 





N(A) / event 


pr slope [GeV l ] 


Data 

CDM 

with P s /P u = 0.3 
with Ps/Pu = 0.2 

MEPS 

with Ps/Pu = 0.3 


0.038 ± 0.006 ± 0.002 

0.066 ± 0.003 
0.050 ± 0.002 

0.068 ± 0.003 


-1.4 ± 0.3 ± 0.1 

-1.04 ± 0.07 
-1.00 ± 0.06 

-0.98 ± 0.06 



Table 3: Results of the A measurement for events with 10 < Q 2 < 640 GeV 2 , 0.0003 < x < 0.01, 
y > 0.04 and for a A with 0.5 < pr < 3.5 GeV, \rj\ < 1.3. The two CDM samples have been generated 
with a different strange-quark suppression factor P s /P u . 
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We have studied the mean K° multiplicity and the ratio of K° to charged particle multiplicities 
N(if°)/N (tracks) as a function of the Q 2 of the event. In order to stay in the region of uniform 
acceptance given by the inner tracking detector geometry and the analysis cuts, we restrict this study 
to events with —1.5 < r/ 7 , t < 0. Figure [| shows the distribution of our event sample in the (x,Q 2 ) 
plane. The lines of constant 7/j delimiting the accepted events and the Q 2 bins chosen for this study 
are shown. In those bins the variables Q 2 and W are correlated: as Q 2 increases from 10 GeV 2 to 
200 GeV 2 , the mean value of W increases from 110 GeV to 160 GeV. Figure ||a, b show the mean 
K° multiplicity and the ratio N(i*r )/N(tracks) in the selected bins plotted versus the mean Q 2 of the 
bins. The number of charged particles does not include secondary particles from K° and A decays 
and from weakly decaying particles with a lifetime > 10~ 8 s. A slight increase of the K° multiplicity 
and a constant behaviour of N(isT°)/N(tracks) are observed. We have included the predictions from 
the CDM and MEPS Monte Carlo samples, which describe the data reasonably well. A study at 
the Monte Carlo generator level shows that the mean K° multiplicity is independent of Q 2 for fixed 
W . Since data and Monte Carlo agree over a wide range of Q 2 , we conclude that the mean K° 
multiplicity of our data also shows no Q 2 dependence at fixed W within the accuracy of these data. 
Furthermore the ratio of K° to charged particle multiplicities is observed to be constant and thus 
within our experimental errors this ratio does not depend on the kinematic variables in the region 
under study. Therefore we attribute our observed increase of K° multiplicity with Q 2 to the increase 
of the corresponding W values. 

8.2 K° production in events with a large rapidity gap 

The DIS data sample is a mixture of non-diffractive and diffractive events. We have searched for 
differences in K° production in these event types. Following our earlier publications [36^ ], we 
separate a non-rapidity gap event sample (NRG) and a LRG event sample using ?? max . 7? max is the 
largest pseudorapidity of any calorimeter cluster in an event, where a cluster is defined as an isolated 
set of adjacent cells with summed energy above 400 MeV. The NRG sample is selected by ry max >1.5. 
It is dominated by non-diffractive events. The requirement rj mauX < 1.5 selects a LRG sample which 
is dominated by diffractive events. The standard non-diffractive DIS models (CDM, MEPS) give a 
reasonable description of the r/ max distribution for values above 1.5 but cannot account for the excess of 
events at lower values (see Fig. 0a). Values of ry max > 4.3, which are outside the calorimeter acceptance, 
occur when energy is deposited in many contiguous cells around the beam pipe in the proton direction. 
An admixture of about 10% - 20% of diffractive events generated with the NZ or POMPYT Monte 
Carlo programs to the non-diffractive Monte Carlo sample gives a reasonable description of the /j max 
distribution. The background of non-diffractive DIS events in the LRG sample is estimated to be 7% 
||. Less than 10% of the NRG DIS event sample are diffractive events. Figure ^b shows the rj max 
distribution for those events which have a candidate in the signal band. The r? max distribution 
of events from one of the non-diffractive (CDM) and from one of the diffractive (NZ) Monte Carlo 
samples is also shown. The excess of Kg candidates over predictions from the CDM model for ?7ma X 
< 1.5 represents the K s production in diffractive events. 

As discussed elsewhere || |3(| , the acceptances for diffractive events selected by the LRG requirement 
(j/ max ) and for NRG events are flat with respect to W and Q 2 for W > 140 GeV. We have therefore 
restricted our comparison to events with W > 140 GeV. After this additional requirement, 11000 NRG 
events and 940 LRG events remain. In the LRG sample we find in the signal region 18 K s candidates 
over a background of 2 candidates. 

Figure || shows the differential K° multiplicity as a function of the transverse momentum and pseudo- 
rapidity for NRG and for LRG DIS events separately. The results in this subsection are not corrected 
for either the r/ max or the W selection criteria. The predictions of the standard DIS Monte Carlo 
programs (CDM and MEPS) and the diffractive DIS Monte Carlo programs (POMPYT and NZ) are 
shown. The pt distributions have similar shapes in both data subsamples, although the multiplicity 
is lower for the LRG DIS events. Within the limited statistics of the data, both diffractive models 
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give a reasonable description of the K° multiplicities in LRG events. 

Since the invariant mass of the measured hadronic system in LRG events is smaller than in NRG 
events, a reduced K° rate is expected in the diffractive events. We have compared the K° multiplicity 
with the charged particle multiplicity for both subsamples. Table |] lists the K° multiplicity and the 
ratio of the K° to charged particle multiplicity for NRG and LRG DIS events and for the Monte Carlo 
samples. If one subtracts the diffractive background, which, as seen from Fig. ^, is still present in the 
NRG DIS sample, the quoted K° multiplicity in the non-diffractive DIS sample increases by 5%. The 
ratios of K 0, s to charged tracks for both data samples are consistent with each other. Thus, within 
the limited statistics, these results give no indication of any additional strange quark enhancement or 
suppression in the production mechanism of the LRG final state. 





Data type 


N(i^ u ) /event 


N(K U ) / N(tracks) 


NRG 


ZEUS data 

CDM 

with P s /P u = 0.3 
with P s /P u = 0.2 

MEPS 

with P s /P u = 0.3 


0.344 ± 0.023 ± 0.025 

0.396 ± 0.009 
0.296 ±0.011 

0.375 ± 0.009 


0.077 ± 0.006 ± 0.008 

0.095 ± 0.003 
0.071 ± 0.003 

0.096 ± 0.003 


f/max <l-5 

LRG 


ZEUS data 

POMPYT 

NZ 


0.156 ±0.047 ±0.007 
0.106 ±0.010 
0.173 ±0.017 


0.071 ±0.021 ±0.007 
0.058 ± 0.006 
0.073 ± 0.007 



Table 4: The K° multiplicity and the ratio of the K° and charged particle multiplicities for NRG 
and LRG DIS events. The predictions of five Monte Carlo samples are listed. The diffractive samples 
are generated with a strangeness suppression factor P s /P u = 0.3. 



9 Study of systematic errors 

We have investigated several sources of systematic errors for our measurements of the K° and A 
production rates. 

1) The sensitivity of the results with respect to the track and primary vertex reconstruction methods 
was determined by repeating the analysis with a modified version of the reconstruction package. The 
differences seen are at the 5% level for the multiplicity distributions. No systematic effect is apparent. 
The ratio of K° to charged particle multiplicity is similarly unaffected. 

2) The sensitivity of the results on the choice of the Kg and A selection criteria has been investigated 
by varying them by ±25% of their nominal values. The uncertainty in the results from the DIS event 
selection was determined by repeating the analysis with different electron finding algorithms and by 
varying the event selection criteria by reasonable values. The systematic error from those sources is 
about 5% except for the highest r] and pi points in the multiplicity distributions and for the results 
of the LRG event analysis, where the error is up to 15%. The mean particle multiplicities per event 
show lower systematic errors (3%) than the bin by bin errors in the figures. 

3) Uncertainties from events rejected by the DIS event selection criteria and event migration effects 
were determined by detailed Monte Carlo studies of the K° and A production in the events migrating 
into and out of the selected Q 2 ,x,y range. The K° and A rate of events migrating into this range is 
comparable to that of events migrating out. The uncertainty from these sources is at the 5% level. 
The additional kinematic restriction of W > 140 GeV for the LRG comparison introduces a higher 
uncertainty (7%) for the results. The mean particle multiplicities show a 2% uncertainty for NRG 
DIS events and 5% for the LRG DIS events. 
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4) We determined a photoproduction contamination in the event sample of 2.5%. The event sample 
which was kinematically restricted to W > 140 GeV contains a higher background of 3.5% as shown 
in |L8], . We have estimated how these photoproduction events affect our analysis by studying the 
stability of the results when varying the scattered electron energy and the 5 selection criterion. We 
quote an uncertainty from this source of 3%. The influence on the results from initial state radiative 
events not removed by the 5 requirement is below 3% except for the lowest r\ point in Fig. ||b where 
it is 15%. 

5) The K° multiplicity versus Q 2 is rather sensitive to the background below the M nn signal. The 
combinatorial background increases with Q 2 due to the observed higher particle multiplicity in events 
with higher Q 2 . Also migration effects are non- negligible. Both effects together may induce variations 
to the measured values between -11% and +3% depending on the Q 2 bin and on the Monte Carlo 
simulations used to determine them. We include an overall systematic error of 10% to our results from 
these sources. 

6) The results for the ratio of the K° multiplicity to the charged particle multiplicity are affected by 
uncertainties similar to those for the K° multiplicity alone. The variations resulting from different 
correction procedures of calculating the mean charged track multiplicity or from using different Monte 
Carlo samples for the correction are within a few percent. The relative changes of the ratios of K° to 
charged particle multiplicities for the NRG and the LRG data samples are below 5% when different 

ranges for the charged particles are considered. 

7) The strange quark density of the proton structure function does not affect our acceptance correc- 
tions. 

10 Summary and discussion 

We have measured the K° and A multiplicities for deep inelastic ep scattering events at yfs = 296 GeV 
with 10 GeV 2 < Q 2 < 640 GeV 2 , 0.0003 < x < 0.01 and y > 0.04 in the ZEUS experiment at 
HERA. We have restricted the analysis to the K° and A kinematic region p^ > 0.5 GeV and < 1.3. 
About 23% (20%) of the K° (A) are predicted to be produced within this kinematic range. 
In this kinematic range the mean number of K° (A) per event is 0.289 ± 0.015 ± 0.014 (0.038 ± 0.006 
± 0.002). The results on particle production from lower energy e + e~ data, which are incorporated in 
the current DIS Monte Carlo simulation programs (i.e., strange quark suppression factor P s /P u = 0.3), 
predict higher K° and A multiplicities than those observed in the data. Using a smaller value of 
0.2 reduces the predicted multiplicity and gives a better agreement with the data, especially for A 
production. Nevertheless, with P s /P u = 0.2 the prediction for the A multiplicity is still higher, while 
the prediction for the K° multiplicity is lower than the measured values. The Monte Carlo models allow 
an adjustment of the production rates of the different particle types by changing other parameters, like 
the ratio of diquarks to single quarks created from the sea, P qq /P q , as well as the suppression factor 
for strange diquarks, (P us / P u d) I (Ps / Pd) ■ Our results indicate the need for tuning these parameters 
which requires a detailed measurement of the ratios of pions, kaons, lambdas and protons over a larger 
kinematic range. This is beyond the scope of this paper. The shapes of the distributions for K 0, s and 
A's are described by both models and do not depend on the chosen parameter P s /P u . 
The mean K° multiplicity of our data shows no indication for a Q 2 dependence at fixed W. Also, the 
ratio of K to charged particles is observed to be independent of the kinematic variables in the range 
studied. 

We observe K° production in DIS events with a large rapidity gap with respect to the proton direction. 
The K° multiplicity in LRG events is approximately a factor of two lower than in non-diffractive DIS 
events. The ratio of K° to charged particles is found to be the same in both samples. Thus we observe 
no additional enhancement or suppression of neutral kaon production in events with a large rapidity 
gap compared to events without a gap. 
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Figure 1: Invariant masses for two particle combinations assuming a pn versus a tttt mass hypothesis. 
The background to the Kg signal is reduced by removing the candidates with M p7r < 1.12 GeV. 
The background to the A signal is reduced by removing candidates with M nn > 0.481 GeV. The 
dashed lines correspond to these mass values. Candidates for photon conversions removed by requiring 
M ee < 50 MeV are indicated by squares. 
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ZEUS 1993 




Figure 2: a) Reconstructed K s mass distribution and b) reconstructed A mass distribution in accepted 
events with 10 < Q 2 < 640 GeV 2 , 0.0003 < x < 0.01, and y > 0.04. The kinematic range of accepted 
Kg candidates is 0.5 < pr < 4 GeV and for A candidates it is 0.5 < pr < 3.5 GeV. The rj range of 
accepted Kg and A candidates is \rj\ < 1.3. 
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Figure 3: a) Differential multiplicity of K° versus transverse momentum of the kaons in the restricted 
kinematic rj and range; b) same for the pseudorapidity of the kaons. The inner error bars show 
the statistical error and the outer ones correspond to the statistical and systematic errors added in 
quadrature. The predictions of the three models discussed in the text are shown. 
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Figure 4: a) Differential multiplicity of A versus transverse momentum in the restricted kinematic pr 
and v) range; b) same for the pseudorapidity of the A's. The inner error bars show the statistical error 
and the outer ones correspond to the statistical and systematic errors added in quadrature. Statistical 
errors dominate. The predictions of the three models discussed in the text are shown. 
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DA 

Figure 5: The (x, Q 2 ) plane with lines of constant r\ lh . For the determination of the Q 2 dependence 
of the K° multiplicity only events with -1.5 < r/ 7h < are accepted. The horizontal lines correspond 
to the Q 2 bins chosen. 
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Figure 6: a) Mean multiplicity of K° versus the Q 2 of the event in the restricted kinematic ranges, b) 
Ratio of the mean K° multiplicity to the mean charged particle multiplicity as a function of the event's 
Q 2 . The inner error bars show the statistical error and the outer ones correspond to the statistical 
and systematic errors added in quadrature. The predictions of three Monte Carlo samples are shown. 
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Figure 7: a) r/ max distribution for all selected events, for one of the non-diffractive Monte Carlo 
samples (CDM) and for one of the diffractive Monte Carlo samples (NZ). The relative fractions of 
CDM (88%) and NZ(12%) events are chosen so that their sum reproduces the data distribution best, 
b) ??max distribution of events with a K s candidate in the signal region of the M W7T distribution. The 
predictions of the CDM model and of the NZ model are overlaid. The dotted line corresponds to r] max 
= 1.5. 
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Figure 8: Differential K° multiplicity in NRG DIS and LRG DIS events as a function of their a) 
transverse momentum; and b) pseudorapidity. The inner error bars show the statistical error and the 
outer ones correspond to the statistical and systematic errors added in quadrature. The predictions 
of five Monte Carlo samples are shown. The two diffractive Monte Carlo samples are generated with 
Ps/Pu = 0.3. 
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